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Abstract 

It is currently possible to fabricate crystalline silicon solar cells with the absorber thickness 
ranging from a few hundreds of micrometers (conventional wafer-based cells) to devices as thin 
as 1 /mi. In this work, we use a model single-junction solar cell to calculate the limits of energy 
conversion efficiency and estimate the optimal absorber thickness. The limiting efficiency for cells in 
the thickness range between 40 and 500 pin is very similar and close to 29%. In this regard, we argue 
that decreasing the thickness below around 40 pm is counter-productive, as it significantly reduces 
the maximum achievable efficiency, even when optimal light trapping is implemented. We analyse 
the roles of incomplete light trapping and extrinsic (bulk and surface) recombination mechanisms. 
For a reasonably high material quality, consistent with present-day fabrication techniques, the 
optimal thickness is always higher than a few tens of micrometers. We identify incomplete light 
trapping and parasitic losses as a major roadblock in improving the efficiency upon the current 
record of 25.6% for silicon solar cells. Finally, considering the main parameters that impact solar 
cell performance, we quantify the constraints and requirements for achieving a specified energy 
conversion efficiency, which is important for a proper design strategy of high efficiency silicon solar 
cells. 


* piotr.kowalczewski@unipv.it 


1 



I. INTRODUCTION 


In this work we focus on the efficiency limits of crystalline silicon (c-Si) solar cells. In 
this regard, we consider both ideal devices (perfect material and interfaces), as well as 
more realistic conditions, including defect-related recombinations, parasitic losses, and non- 
optimal light management. It is currently possible to fabricate c-Si solar cells with absorber 
thickness values that differ by two orders of magnitude. On the one end of the thickness 
range there are conventional wafer-based c-Si solar cells with the absorber thickness of the 


order of a few hundreds of micrometers jl]. Epitaxial growth allows fabricating solar cells 
with the thickness of a few tens of micrometers Q. Finally, epitaxy-free fabrication { 3 , 4\ 
makes possible fabricating c-Si solar cells with the absorbing layer as thin as 1 /mi. For a 
proper design strategy, it is important to determine the optimal absorber thickness in terms 
of the energy conversion efficiency. We address this question using a model single-junction 
silicon solar cell. 


It is well know that reducing the silicon thickness leads to the reduction of total ab¬ 
sorption, and thus of the photocurrent. This has to be compensated by implementing an 
appropriate light-trapping scheme. Yet, it should be emphasized that even when the optimal 
light trapping is applied (i.e., perfect anti-reflection action combined with a Lambertian scat- 
terer), the maximum achievable absorption still decreases with decreasing material thickness 
( 5 , 6]. On the other hand, V oc generally tends to decrease with increasing thickness. For 
a given material quality, this leads to an optimal thickness that maximizes the conversion 
efficiency js, 7]. In this context, the goal of this work is twofold: First, to establish the 
most suitable thickness range for c-Si solar cells to approach the efficiency limits. Second, to 
identify the parameters that have to be improved in order to increase the energy conversion 
efficiency beyond the current record of 25.6% [§]. 


We use efficiency as a figure of merit to assess different solar cell structures. This is 
motivated by the fact that the cost of electricity is mainly determined by the efficiency 
rather than by the cost of the active material (which, in the case of silicon, is constantly 
decreasing). I 11 this regard, it is particularly important to estimate the optimal absorber 
thickness range that maximizes the efficiency. In our analysis we include intrinsic Auger 
recombination, as well as defect-based bulk and surface recombination mechanisms. We 
also introduce a simple approach that allows us to consider parasitic losses. We compare 
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our results with the measured performance of state-of-the-art silicon solar cells, showing 
the room for improvement in terms of current, voltage, and fill factor. We pay particular 
attention to parasitic losses, which we believe are a major roadblock in reaching efficiency 


above the current record of 25 


8 ], 


A review of the efficiency limits of silicon solar cells is given in Ref. [9]. The limits 


& 


reported in the literature 10-|l2] are usually calculated using the idealized diode equation 


13]. This approach has the following limitations: 


1. The diode equation gives less accurate results when the cell thickness is decreasing. We 
attribute this inaccuracy to the assumptions underlying the treatment of the space- 
charge region (SCR) in the idealized diode formalism. In practice, idealized diode 
equation tends to significantly overestimate efficiency for thin cells. We further discuss 
the accuracy of the results obtained using the idealized diode equation in the Appendix. 

2. Solar cells require selective contacts, which can be achieved using a p-n junction. Yet, 
the junction is not explicitly considered in the ideal diode approach. Therefore, ideal 
diode equation gives unrealistic results that overestimate efficiency for undoped silicon. 


In our approach, we calculate the efficiency limits overcoming these limitations: we analyt¬ 
ically obtain the photogeneration rate (ranging from double-pass absorption to Lambertian 
light trapping), explicitly consider a p-n junction, and numerically solve the drift-diffusion 
equations. This allows us to calculate more realistic efficiency limits of silicon solar cells in 
a wide range of cell thicknesses and doping levels. 

The paper is organized as follows: in Sec. HT]we describe our approach, based on analytical 
photogeneration rate and numerical solution of drift-diffusion equations. In Sec. IIIII we cal¬ 
culate the efficiency limits for c-Si solar cells. In Sec. IIVI we discuss the effects of incomplete 
light trapping in approaching the efficiency limits. In Sec. [V] we cover the impact of bulk 
material imperfections on the cell performance, and introduce a simple approach to include 
parasitics losses in the analysis. In Sec. ED we discuss the role of surface recombination. 
In Sec. IVIII we quantify the requirements, in terms of bulk and surface material quality, to 
achieve a given efficiency level. Conclusions are given in Sec. IVIII! Finally, in the Appendix 
we compare the results of our numerical treatment with those obtained using the idealized 
diode equation. 
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FIG. 1. Investigated solar cell structure consisting of a 5nm thick n-type emitter (the blue region), 
p-type base (the green region), and a perfect back reflector (BR), which also serves as a back 
contact. The front surface is textured so that the incident light is scattered and trapped within 
the absorber. 

II. NUMERICAL APPROACH 


Let us consider the structure sketched in Fig. |T[ It consists of a 5 nm thick n-type emitter 
and a p-type base of variable thickness. Such a thin emitter allows to minimise recombination 
losses in this heavily doped layer. To calculate the efficiency limits we assume no reflection 
at the front interface, a perfect back reflector (BR), and a Lambertian light trapping QQ- 
In Fig. |T]we schematically show the Lambertian scatterer at the front, yet we note that the 
photogeneration rate is calculated analytically and in the electrical calculations the structure 
is assumed to be flat. Finally, we assume full-area contacts: the carriers are collected at the 
silicon/BR and emitter/front contact interfaces. 

The photogeneration rate corresponding to the Lambertian limit is calculated as in 
Ref. |(3|: 


^ ^ a lt (R h e~ 2 ^ w e^ z + e~^ z ) 

GlL{ZiE) — - 7 -—- X 0AM1.5G) 


1 - e~ 2oiltW 1 


responding to the AM1.5G solar spectrum 


( 1 ) 


where nsi is the refractive index of c-Si [16] and 0 ami.5 G is the photon flux density cor- 

17]. The effective absorption coefficient cp t is 
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related to the light path enhancement in textured cells according to Ref. 15]. In the liter¬ 


ature we can find a number o: 
limit, including ordered 


random textures 


24 


18 


different strategies that allow approaching the Lambertian 


20] and quasi-ordered 21-]23] photonic structures, as well as 


;2 71 . 


The photogeneration rate calculated using Eq. |Tj is integrated with respect to energy 
(over the solar spectrum) and used as the generation term in the drift-diffusion equations, 
which are solved using Finite-Element Method FEM). We use FEM implementation in 


the commercial device simulator Si 


described in our previous works 


t 


vaco Atlas 


28]. This methodology is similar to the one 


29]. 


III. EFFICIENCY LIMITS 


We start by considering solar cells limited by intrinsic Auger recombination. We treat 
Auger recombination using the parametrization reported in Ref. 30] - We also consider 
band gap narrowing (BGN) according to the model by Schenk 31] - Finally, we neglect free 
carrier absorption, which is a second-order effect [101, 12]. Silicon is an indirect band gap 
material, and thus we also neglect losses related to radiative recombination, which may give 
an appreciable effect only for very thick cells. Yet, for the thick cells the probability of photon 
recycling increases [ill, i.e., radiatively emitted photons are reabsorbed. These two effects 
arc likely to compensate each other. In the appendix we show that radiative recombination 
together with photon recycling have a negligible effect on the cell performance. 

In Fig. [2](a) we show the limiting efficiency of c-Si solar cells as a function of the absorber 
thickness, calculated for the structure sketched in Fig. [H For each thickness, we have 
simultaneously optimized the emitter and base doping concentrations (the doping profile in 
each layer is assumed to be constant). The efficiency as a function of emitter doping N ,d has 
a wide maximum, and the optimal value N ,d = 1.5 x 10 18 cm -3 does not change with the 
absorber thickness, as the emitter thickness itself is kept constant. Regarding the base doping 
N a , the solid lines in Fig. [2] are calculated by optimizing N a for each thickness. On the other 
hand, the dashed lines are calculated assuming a constant doping N a = 10 16 cm^ 3 , which is 
the optimal doping for the optimal absorber thickness equal to 170 pm. The optimal base 
doping is highest for thinner cells (reaching 10 17 cm -3 for the 1 pm thick cell) and decreases 
with increasing thickness. 
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FIG. 2. (a) The limiting efficiency of c-Si solar cells as a function of the absorber thickness, (b) 
V oc and J sc , and (c) fill factor (FF) corresponding to the calculated efficiency limits. The solid 
lines are calculated by optimizing the base doping for each thickness, whereas the dashed lines are 
calculated assuming the constant doping N a = 10 16 cm~ 3 . The triangles denote the performance 
of the record-efficiency HIT cells [32]. 
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FIG. 3. Efficiency as a function of the base doping and absorber thickness. The dashed line 
indicates the optimal doping. 


The maximum efficiency is equal to ?/ max = 28.73%, and it is obtained for the 170 /irn 
thick cell. This is lower than the efficiency limit of 29.43% reported recently in the literature 


12]. The main reason is that in our calculations, the base doping increases Auger losses and 


therefore reduces the efficiency. In the contrary, the limit reported in the literature has been 
obtained using the idealized diode equation and assuming undoped silicon, as explained in 
Sec. HI and in the Appendix. 

We note that the efficiency as a function of the absorber thickness has a broad maximum, 
and even a small change in the input parameters can substantially shift the nominal optimal 
thickness. In the range between 40 and 500 /irn, the calculated efficiency differs by no more 
than 1% (relative units) from r/ max . This somehow arbitrary interval shows that for thick¬ 
ness values that differ considerably (by one order of magnitude), the maximum achievable 
efficiency is very similar. Finally, for an absorber thickness below 40 /irn, the efficiency drops 
significantly. 

These results are compared with the efficiencies of the state-of-the-art HIT cells 32] (the 
triangles). We note that the efficiency of the HIT solar cell has been recently increased to 


B- 


25.6 % |8fl. Yet, the thickness of the record cell is not specified, and therefore we were unable 
to include it in Fig. [2j This analysis shows that the theoretical margin for improvement of 
single-junction silicon solar cells is around 3% (absolute units). 
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In Fig. |2fb) we show V oc and J sc corresponding to the calculated efficiency limits. The 
base doping optimization allows to slightly improve V%, whereas J sc is practically the same. 
Comparing these results with the performance of the HIT cells, we can see that for both 
values of thickness (98 and 151 pm), the measured value is around 99% of the limiting V oc , 
around 94.6% of the limiting FF, and around 91.6% of the limiting J sc . Therefore, these 
cells are limited by, in decreasing order of significance, J sc , FF, and V oc . This shows that 
optimizing light trapping and, especially, minimizing optical parasitic losses (affecting J sc ), 
as well as improving the contacts (affecting FF) is the key to further increase the efficiency 
towards the limiting values. 

To further elaborate the issue of the base doping optimization, in Fig. [3] we show the 
efficiency as a function of the base doping and absorber thickness. This plot can be com¬ 
pared with the results presented in Fig. 4 reported in Ref. [ijj], which are obtained using 
the idealized diode equation. In the case of the diode equation, solar cells approach the 
limiting efficiency for a practically undoped silicon. Yet, in our calculations we can see a 
clear maximum around 10 16 cm -3 (for the thickness range around 100 - 200pm). Then, 
efficiency decreases with decreasing doping. This difference is due to the fact that the junc¬ 
tion is not explicitly considered in the ideal diode approach. Therefore, the requirement of 
selective contacts, which can be achieved using a p-n junction, is not included. This leads 
to unrealistic results for very lightly doped materials. 

For the simplicity of the analysis, in the rest of this paper we assume the constant base 
doping !V a = 10 16 cm~ 3 . 


IV. EFFECTS OF INCOMPLETE LIGHT TRAPPING 

In the calculations above we have assumed a Lambertian light trapping, which is often 
taken as a benchmark in the optical design of solar cells. Yet, it is difficult to fulfil this 
assumption in realistic devices. Therefore, let us now focus on the role of light trapping 
in achieving the efficiency limits. To do so, we consider solar cells with incomplete light 
trapping: the photogeneration rate is taken as a weighted average of the photogeneration G* 2 P 
corresponding to the double-pass absorption (i.e., unstructured cell) and the photogeneration 
G*ll corresponding to the Lambertian limit. G 2p is calculated as 
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( 2 ) 


G 2p (z, E) = a (e~ za + x 0 AM i.5G, 


where a is the absorption coefficient of c-Si 16]. As previously, also in the double-pass 
case we assume a perfect back reflector and anti-reflection action. We note that there 
are no interference effects included in Eq. ([2]). When we consider single wavelengths, the 
interference effects are profound. Yet, when we integrate the photogeneration rate over solar 
spectrum, the interference peaks are smeared out, which justifies the approximation used in 
Eq. (J2D- The resulting weighted photogeneration rate is calculated as 


G(z) = P\ oss X [(1 - LF) x G 2p (z) + LF x G LL (z )}, 


(3) 


where LF is the light-trapping factor: LF = 0 corresponds to the double-pass case, whereas 
LF = 1 corresponds to the Lambertian limit. G is a function of depth z. Finally, Pi oss 
is a factor related to reflection at the front interface and to parasitic optical losses. For 
simplicity, we assume that Pi oss does not depend on wavelength. 

We begin by assuming no parasitic losses, that is Pi oss = 1. In Fig. [Ha) we demonstrate 
that light trapping increases the maximum achievable efficiency. This is not the case for 
semiconductors with a direct energy band-gap, like GaAs, where surface texturing does not 
increase the maximum efficiency 33]. Therefore, for silicon solar cells light trapping is an 
essential element required to approach the efficiency limits, even for very thick cells. 

As expected, light trapping significantly increases J sc , which is shown in Fig. EOT). J sc 
for the thickest cells does not saturate because of the small absorption below the energy 
band-gap. Moreover, Fig. [4](c) shows that light trapping also slightly improves V oc . Yet, this 
effect is appreciable only for cell thicknesses below 10 /mi, that is in a thickness range which 
is not particularly promising for achieving high efficiency. V oc as a function of the thickness 
for the double-pass case exhibits a gentle maximum. This may be partly because the base 
doping is not optimized for thin cells. 

Although the importance of light trapping for maximizing the conversion efficiency of 
c-Si solar cells is well known, the present results allow to quantify the effects of incomplete 
light trapping on the efficiency. Moreover, Eq. (J3l) represents a simple model that can be 
also used to study the effects of parasitic losses, which are different from the effects related 
to incomplete light trapping. We shall discuss it in more details in the next section. 
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FIG. 4. (a) Efficiency, (b) short-circuit current J sc , (c) open-circuit voltage V oc calculated as a 
function of the absorber thickness for different values of the light-trapping factor LF. 

V. MATERIAL IMPERFECTIONS AND PARASITIC LOSSES 

So far, we have considered an idealized material, and therefore the cell performance was 
limited by intrinsic Auger recombination. Let us now consider extrinsic losses related to 
defect-base SRH recombination and parasitic optical losses. 

With the emitter as thin as 5 nm, the cell efficiency is likely to be limited by the diffusion 
length L n ,SRH of the minority carriers (electrons) in the base. In Fig. [5](a) we show the 
efficiency as a function of the absorber thickness and L nj srh, calculated for the structure 
with Lambertian light trapping. The diffusion length of holes in the highly-doped emitter is 


10 













taken equal to L Pi srh = L n / 10. It can be seen that thicker cells are more sensitive to SRH 
recombination. In this regard, thinner cells with the absorber thickness of around 40 /irn 
can reach nearly the same ?/ max as the thicker cells, but the conditions in terms of absorber 
quality are relaxed. Finally, reducing the absorber quality (i.e., decreasing L^srh) shifts the 
optimal thickness towards thinner cells. 

Results shown in Fig. [5|b) are calculated including parasitic losses Pi oss = 0.9. We 
estimate that this is approximately the level of parasitic losses in the HIT cells discussed 
above. The trends are similar to those calculated assuming P\ oss = 1, yet the maximum 
achievable efficiency is proportionally decreased. Moreover, for a given efficiency level, the 
optimal thickness is increased. 

The results of Fig. fallow us to estimate the optimum thickness and the required material 
quality to reach a given efficiency level. For example, to surpass 26% efficiency, in the 
idealized case of no parasitic losses and complete (Lambertian) light trapping, we need an 
electron diffusion length around 1 mm. In the more realistic case, when parasitic losses are 
included (say, Pi oss = 0.9), achieving efficiency above 26% requires a diffusion length higher 
than fa 4 mm for an optimal thickness of 50 pm. 
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FIG. 5. (a) Efficiency as a function of the thickness and diffusion length L nj sRH of the minority 
carriers (electrons) in the base, calculated for the structures with a Lambertian light trapping, (b) 
The same quantity but including factor related to parasitic losses Pi oss = 0.9. We note that L n ,SRH 
is related to SRH recombination, originating from defects in the base. The scale on the right shows 
the corresponding carrier lifetime, calculated assuming electron mobility ji n = 1548 cm 2 V -1 s -1 . 
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FIG. 6. Efficiency as a function of (a) front and (b) back surface recombination velocity (SRV), 
and of the absorber thickness. For the bulk transport losses, we assume only Auger recombination. 

VI. SURFACE RECOMBINATION 

Let us now estimate the constraints on the efficiency imposed by surface recombination, 
which is another extrinsic loss mechanism related to defect states at the surface. In Fig. [H] 
we show the efficiency as a function of (a) front and (b) back surface recombination velocity 
(SRV), and of the absorber thickness. At first, for the bulk transport losses we assume only 
Auger recombination. It can be seen that approaching the efficiency limits requires SRV to 
be less than a few cm/s, with solar cells being more sensitive to recombination at the rear 
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FIG. 7. Efficiency as a function of the front surface recombination velocity (SRV) and the absorber 
thickness. For the bulk transport losses, we assume Auger and SRH recombination. The material 
quality (i.e., L n ,SRH) is reduced, so that the maximum achievable efficiency is 27% for L n ~ 2200 pm 
(a) and 25% for L n ss 840 prn (b). 

interface rather than to recombination at the front. 

In the calculations above, we assumed a perfect material quality. Yet, in our previous 
work we have demonstrated that the importance of surface recombination depends on bulk 
recombination rate, i.e., the higher is the material quality, the more important are losses at 
the surface Q. From now on we focus on front surface recombination, neglecting recom¬ 
bination at the back, which has anyway to be of the order of a few cm/s to target high 
efficiency. In realistic devices, there is usually an oxide passivating layer at the back and 
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point contacts are implemented. For this reason, back SRV can be significantly reduced. On 
the other hand, front SRV may be increased due to the texturing. 

In Fig. |7] we plot the efficiency as a function of the front surface recombination velocity 
(SRV) and absorber thickness. For bulk transport losses, we assume Auger and SRH recom¬ 
bination. We adjust the diffusion length L n g rh of the minority carriers (electrons) in the 
base, so that the maximum achievable efficiency is reduced to a certain value. This kind of 
analysis allows to answer the following question: Having a given material quality, what is 
the maximum allowed value of the front SRV to achieve the limiting efficiency? 

In Fig. [3(a) the maximum possible efficiency is around 27%, which is obtained for L n ps 
2200 /irn. With this material quality, the optimal absorber thickness is 40 /im, and the 
maximum efficiency can be achieved for the thickness range between 20 and 80 /irn. In this 
case, front SRV should be below a few tens of cm/s. 

In Fig. [3(b) we further reduce diffusion length to L n ps 840 /irn, so that the maximum 
possible efficiency is 25%. The region of highest efficiency (25% or slightly above) is clearly 
wider than in the previous plot. Also the conditions for surface recombination at the front 
interface are relaxed: in this case, front SRV should be below a few hundreds of cm/s. 
This analysis confirms that the impact of surface recombination decreases with decreasing 
material quality. 

VII. REQUIREMENTS FOR APPROACHING THE EFFICIENCY LIMITS 

To summarize the results presented in this paper, in this section we estimate the require¬ 
ments for approaching a given energy conversion efficiency. In this regard, we analyse both 
bulk material quality (related to diffusion length of electrons in the base T^srh) and quality 
of the surfaces (related to surface recombination velocity). 

In Fig. |8](a) we show the diffusion length of electrons in the base T n ,SRH as a function of 
the maximum achievable efficiency, whereas Fig. [3(b) indicates the corresponding (optimal) 
thickness of the absorbing layer. These results are extracted from Fig. [3 i.e., for each T n ,SRH 
we have extracted the maximum achievable efficiency and the corresponding cell thickness. 

In Fig. [3(a) we can distinguish two trends. In the case of no parasitic losses (black line), 
for the efficiency range below around 28%, L n and the corresponding maximum achievable 
efficiency scale linearly. Yet, above 28%, significant increase of L n sRH gives only a minor, if 
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FIG. 8. (a) Diffusion length of electrons in the base L n as a function of the maximum achievable 
efficiency, (b) Corresponding (optimal) thickness of the absorbing layer. 


any, improvement of efficiency. 

The red line in Fig. Efa) shows the results calculated including parasitic losses, that is 
floss — 0.9. Parasitic losses result in a shift of the red curve towards lower efficiencies, with 
the maximum efficiency of around 26%. Indeed, the measured record efficiency is slightly 
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FIG. 9. (a) Front surface recombination velocity (SRV) as a function of the maximum achievable 
efficiency, (b) Corresponding optimal thickness of the absorbing layer. We consider cells limited 
by Auger recombination (the black lines), as well as cells with SRH recombination: L n ~ 1700 fim 
with maximum achievable efficiency around 27% (the red lines) and L n ~ 660 ^m with maxi¬ 
mum achievable efficiency around 25% (the blue lines). The colour regions indicate 1% (relative) 
tolerance intervals for each material quality. 
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less than 26%, which identifies parasitic losses as a major roadblock for going beyond 26%. 


We can also see that an optimal absorber thickness changes significantly in the considered 
efficiency range, as shown in Fig. [S(b). For small L n , the optimal thickness is of the order 
of a few tens of micrometers. Then, with increasing L n , the optimal thickness increases 
rapidly, reaching around 170 prn for L n of the order of millimetres. 

A similar analysis can be performed regarding surface recombination. In Fig. |9|a) we 
show front surface recombination velocity (SRV) as a function of the maximum achievable 
efficiency, and in Fig. E^b) we show the corresponding optimal thickness of the absorbing 
layer. These results have been extracted from Fig. [6] and [7) 


We consider cells limited by Auger recombination (the black lines), as well as cells with 
SRH recombination: L n 2200 pm with maximum achievable efficiency around 27% (the 
red line) and L n & 840 pm with maximum achievable efficiency around 25% (the blue lines). 
In general, the cells are fairly insensitive to recombination at the front interface if SRV is 
below a few hundreds of cm/s. This is of the order of magnitude of the effective SRV, 
currently measured for nanostructured solar cells 34], Moreover, increasing SRV increases 
the optimal absorber thickness, as thinner cells are more sensitive to surface recombination. 


As pointed out previously, efficiency as a function of the absorber thickness exhibits a 
wide maximum. For this reason, in Fig. |U(b) we include tolerance intervals for each material 
quality. The intervals are indicated by the colour regions, and show thickness range in which 
efficiency changes by no more than 1% relative. For example, if we consider Auger-limited 
cells (the black lines), and assume front SRV equal to 100 cm/s, we can see that the efficiency 
changes by no more than 1% relative in the thickness range from 80 to 500 pm. We note 
that the calculations are performed up to 500 pm (which we consider a practical limit for 
fabrication). 


The results presented above show that the optimal absorber thickness is determined by 
an interplay between bulk and surface losses. Nevertheless, for a reasonably high bulk and 
surface quality (consistent with present-day fabrication techniques), an optimal thickness 
is always above 100 pm. Therefore, the conclusion that decreasing the absorber thickness 
below a few tens of micrometers is counter-productive holds for a wide range of material 
parameters. However, we emphasize that reducing parasitic losses is required in order to 
increase the energy conversion efficiency above 26%. 
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VIII. CONCLUSIONS 


In conclusion, we presented an electro-optical framework to calculate the limiting ef¬ 
ficiency of c-Si solar cells. In our approach, we take the analytical photogeneration rate 
(assuming full or partial light trapping) and numerically solve the drift-diffusion equations 
to obtain the cell performance. This gives a realistic description of carrier dynamics in the 
device and allows easily introducing intrinsic and extrinsic loss mechanisms. Comparison 
with present-day HIT structures suggests that silicon solar cells are limited by (in decreasing 
order of significance) J sc , FF, and V oc . Thus, improving J sc (via light trapping and reducing 
parasitic losses) and FF (by improving the quality of contacts) is a key to further boost the 
performance beyond the current efficiency record. 

In the first part of this contribution, we calculated the efficiency limits of c-Si solar 
cells. Therefore, we focused on the cells limited by intrinsic Auger recombination. We have 
demonstrated that the limiting efficiency as a function of the absorber thickness exhibits 
a wide maximum: 40 pm thick cell can be nearly as efficient as the solar cell with the 
optimal absorber thickness (around 170/mi). It is therefore more practical to consider the 
optimal thickness range, rather than a single optimal thickness. In this regard, we argued 
that decreasing the thickness below around 40 pm is counter-productive, as it significantly 
reduces the maximum achievable efficiency. 


When extrinsic SRH recombination is considered, we notice that thicker cells are more 
sensitive to bulk losses, and therefore the conditions to reach the limiting efficiency for 
thinner cells are relaxed. Yet, including surface recombination shows the opposite trend: 
thicker cells are less sensitive to surface recombination, as the surface-to-volume ratio is 
reduced. Nevertheless, for a reasonably high bulk and surface quality, it remains true that 
decreasing the thickness below a few tens of micrometers is counter-productive in terms of 
the efficiency. 

The current efficiency record for single junction silicon solar cells is 25.6% 8, [35]. Based on 
the analysis presented in this paper, we believe that the efficiency above 26% is a reasonable 
next step, which is within the reach of present-day fabrication techniques. Yet, this requires, 
in the first place, reduction of parasitic losses. Surpassing an efficiency of around 27% puts 
severe constraints on surface recombination. Increasing the efficiency above the calculated 
limit for single junction c-Si solar cells requires novel technologies, like silicon-perovskite 
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tandem structures 
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Appendix A: Accuracy of the idealized diode equation 

The efficiency limits of silicon solar cells are usually calculated using idealized diode 
equation [12]: 


J(V) = J L - qWR, 


(Al) 


where Jl is the photogenerated current density, W is the cell thickness, and R is the total 
recombination rate. In this appendix, we elaborate on the difference between the results 
obtained using Eq. (lAljl and numerical solution of drift-diffusion equations by means of finite- 
element approximation. This allows us to investigate the accuracy of the results obtained 
using Eq. (lAljl in a wide absorber thickness range. 

In Fig. [10] we show (a) efficiency, (b) open-circuit voltage V oc , (c) short-circuit current J sc , 
and (d) fill factor FF as a function of the absorber thickness. The red dashed lines refer to the 
results obtained using the idealized diode equation, whereas the black lines denote the results 
obtained by solving the drift-diffusion equations by means of finite-element approximation. 
The material parameters and recombination models are the same in both approaches: we 
include Auger recombination according to Richter et al. 30] and BGN according to the 
model by Schenk 31]. 

By performing the calculations using Eq. (IA1|) , we aim at reproducing the results reported 
in Ref. 0. In Fig. HOf a) we compare our calculations of efficiency with the results presented 
in the reference work. The results obtained using Eq. (lAljl are very close to the results 


digitalized from Ref. 


12]. We attribute small discrepancies to the fact that in our calculations 


we neglect free-carrier absorption and photon recycling: as discussed above, these two effects 
are likely to compensate each other, and indeed their impact on the cell performance is 
negligible. For this reason, we have also neglected these effects in the calculations presented 
in the main part of this paper. Finally, the calculations presented in this work are performed 
at temperature T = 300 K, whereas the results in the reference work are performed for 
T = 25°C. 

The discrepancy between the results obtained using Eq. (I All) and full numerical sim¬ 
ulations increases with decreasing absorber thickness, as demonstrated in Fig. [lOj In the 
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FIG. 10. (a) Efficiency, (b) open-circuit voltage V oc , (c) short-circuit current J sc , and (d) fill factor 
FF as a function of absorber thickness. The red dashed lines refer to the results obtained using 
the idealized diode equation, whereas the black lines refer to the denote the results obtained by 
solving the drift-diffusion equations by means of finite-element approximation. In (a) we show the 
comparison with the results digitalized from Ref. [12]. 


limiting case of very thick cells, both approaches give essentially the same results. Yet, the 
difference in efficiency for 1 /irn thick cell is close to 2% (absolute value), i.e., the ideal diode 
treatment gives higher efficiency. Since J sc calculated using both approaches is the same [ic]] , 
we can trace the discrepancy back to the difference in V oc and fill factor: V oc calculated us¬ 
ing Eq. (1Alh increases linearly with decreasing thickness, whereas V oc obtained from the full 
numerical simulations tends to saturate. Moreover, FF calculated using Eq. (lAljl increases 
with decreasing absorber thickness, whereas FF obtained from the full numerical simulations 
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shows the opposite trend (small fluctuations may be due to numerical inaccuracy). 

We can therefore conclude that the results obtained using the idealized diode equation 
become less accurate with decreasing absorber thickness, i.e., in this approach the efficiency 
for thin cells is overestimated. We attribute this inaccuracy to the assumptions regarding the 
treatment of the space-charge region (SCR) in the idealized diode formalism: the idealized 
diode equation can be derived from the drift-diffusion equations assuming that there is no 


carrier generation nor recombination in SCR 


41| . It means that the currents in SCR are 


constant. With the assumed doping concentrations, the width of SCR is around 350 nm. 
Therefore, if the cells are thick, SCR is only a small part of the whole cell. Yet, in the case 
of thin cells, SCR is a significant part of the whole cell, and the assumption that there is no 
carrier generation nor recombination in SCR seriously disturbs the current distribution in 
the cell. 


ACKNOWLEDGEMENTS 

This work was supported by the EU through Marie Curie Action FP7-PEOPLE-2010- 
ITN Project No. 264687 ’’PROPHET”. The authors are grateful to Angelo Bozzola for 
carefully reading the manuscript and Marco Liscidini for many fruitful discussions. 


[1] J. Zhao, A. Wang, M. A. Green, and F. Ferrazza, Appl. Phys. Lett. 73 , 1991 (1998). 

[2] J. H. Petermann, D. Zielke, J. Schmidt, F. Haase, E. G. Rojas, and R. Brendel, Prog. 
Photovolt: Res. Appl. 20 , 1 (2012). 

[3] X. Meng, V. Depauw, G. Gornard, O. El Daif, C. Trompoukis, E. Drouard, C. Jamois, A. Fave, 
F. Dross, I. Gordon, and C. Seassal, Opt. Express 20 , A465 (2012). 

[4] C. Trompoukis, O. El Daif, V. Depauw, I. Gordon, and J. Poortmans, Appl. Phys. Lett. 101 , 
103901 (2012). 

[5] A. Bozzola, M. Liscidini, and L. C. Andreani, Optics express 20 , A224 (2012). 

[6] A. Bozzola, P. Kowalczewski, and L. Andreani, J. Appl. Phys. 115 , 094501 (2014). 

[7] P. Kowalczewski, A. Bozzola, M. Liscidini, and L. C. Andreani, J. Appl. Phys. 115 , 194504 
(2014). 


22 







[8] “Panasonic HIT solar cell achieves world’s highest energy conversion efficiency of 25.6% at 
research level,” http://panasonic.co.jp/corp/news/official.data/data.dir/2014/04/ enl40410- 
4/enl40410-4.html (Panasonic Press Release, 10 April 2014, accessed: 28 August 2014). 

[9] R. M. Swanson, in Photovoltaic Specialists Conference, 2005. Conference Record of the Thirty- 
first IEEE (IEEE, 2005) pp. 889-894. 

[10] T. Tiedje, E. Yablonovitch, G. D. Cody, and B. G. Brooks, IEEE Trans. Electron Devices 
31 , 711 (1984). 

[11] M. J. Kerr, A. Cuevas, and P. Campbell, Prog. Photovolt: Res. Appl. 11, 97 (2003). 

[12] A. Richter, M. Hermle, and S. Glunz, IEEE J. Photovolt. 3, 1184 (2013). 

[13] J. Nelson, The Physics of Solar Cells (Imperial College Press, London, 2003). 

[14] E. Yablonovitch, JOS A 72 , 899 (1982). 

[15] M. A. Green, Prog. Photovolt: Res. Appl. 10 , 235 (2002). 

[16] M. A. Green, Sol. Energ. Mat. Sol. Cells 92, 1305 (2008). 

[17] “Reference solar spectral irradiance: Air mass 1.5,” http://rredc.nrel.gOv/solar/spectra/aml.5/. 

[18] A. Abass, K. Q. Le, A. Alu, M. Burgelman, and B. Maes, Physical Review B 85, 115449 

( 2012 ). 

[19] O. Isabella, S. Solntsev, D. Caratelli, and M. Zeman, Prog. Photovoltaics 21, 94 (2013). 

[20] G. Gornard, R. Peretti, E. Drouard, X. Meng, and C. Seassal, Optics Express 21 , A515 
(2013). 

[21] R. Peretti, G. Gornard, L. Lalouat, C. Seassal, and E. Drouard, Physical Review A 88, 053835 
(2013). 

[22] E. R. Martins, J. Li, Y. Liu, V. Depauw, Z. Chen, J. Zhou, and T. F. Krauss, Nature 
Communications 4, 2665 (2013). 

[23] A. Bozzola, M. Liscidini, and L. C. Andreani, Prog. Photovoltaics (2013), 10.1002/pip.2385. 

[24] K. Jager, M. Fischer, R. A. van Swaaij, and M. Zeman, Optics Express 21 , A656 (2013). 

[25] C. Battaglia, M. Boccard, F.-J. Haug, and C. Ballif, Journal of Applied Physics 112 , 094504 

( 2012 ). 

[26] S. Wiesendanger, M. Zilk, T. Pertsch, F. Lederer, and C. Rockstuhl, Applied Physics Letters 
103 , 131115 (2013). 

[27] P. Kowalczewski, M. Liscidini, and L. C. Andreani, Optics Express 21, A808 (2013). 


23 


[28] “Silvaco Atlas Manual,” http://www.silvaco.com/products/ 

tcad/device_simulation/atlas/atlas.html (2014). 

[29] L. C. Andreani, A. Bozzola, P. Kowalczewski, and M. Liscidini, Solar Energy Materials and 
Solar Cells (2014 (in press, doi:10.1016/j.solmat.2014.10.012). 

[30] A. Richter, S. W. Glunz, F. Werner, J. Schmidt, and A. Cuevas, Physical Review B 86, 
165202 (2012). 

[31] A. Schenk, Journal of Applied Physics 84, 3684 (1998). 

[32] M. Taguchi, A. Yano, S. Tohoda, K. Matsuyama, Y. Nakamura, T. Nishiwaki, K. Fujita, and 

E. Maruyama, IEEE J. Photovolt. 4, 96 (2014). 

[33] O. D. Miller, E. Yablonovitch, and S. R. Kurtz, IEEE J. Photovolt. 2, 303 (2012). 

[34] A. Ingenito, O. Isabella, and M. Zeman, Prog. Photovoltaics (2015), 10.1002/pip.2606, 

[35] M. A. Green, K. Emery, Y. Hishikawa, W. Warta, and E. D. Dunlop, Prog. Photovoltaics 
23, 1 (2015). 

[36] O. Malinkiewicz, A. Yella, Y. H. Lee, G. M. Espallargas, M. Graetzel, M. K. Nazeeruddin, 
and H. J. Bolink, Nature Photonics 8, 128 (2014). 

[37] T. P. White, N. N. Lai, and K. R. Catchpole, IEEE Journal of Photovoltaics 4, 208 (2014). 

[38] C. D. Bailie, M. G. Christoforo, J. P. Mailoa, A. R. Bowring, E. L. Unger, W. H. Nguyen, 
J. Burschka, N. Pellet, J. Z. Lee, M. Gratzel, et al., Energy & Environmental Science (2015). 

[39] P. Loper, B. Niesen, S. Moon, S. Martin de Nicolas, J. Holovsky, Z. Rernes, M. Ledinsky, 

F. Haug, J. Yum, S. De Wolf, et al ., (2014). 

[40] The fact that J sc is the same in both cases also confirms that the mesh used in the finite- 
element calculations is dense enough to properly reproduce the photogeneration profile. 

[41] N. W. Ashcroft and N. D. Mermin, Solid State Physics (Brooks/Cole, 1976). 


24 


